Model solution for volume reflection of relativistic particles in a bent crystal 
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For volume reflection process in a bent crystal, exact analytic expressions for positively- and 
negatively-charged particle trajectories are obtained within a model of parabolic continuous potential 
in each interplanar interval, with the neglect of incoherent multiple scattering. In the limit of the 
crystal bending radius greatly exceeding the critical value, asymptotic formulas are obtained for the 
particle mean de?ection angle in units of Lindhards critical angle, and for the final beam profile. 
Volume re?ection of negatively charged particles is shown to contain effects of rainbow scattering 
and orbiting, whereas with positively charged particles none of these effects arise within the given 
model. The model predictions are compared with experimental results and numerical simulations. 
Estimates of the volume re?ection mean angle and the final beam profile robustness under multiple 
scattering are performed. 
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I. INTRODUCTION 

The volume re?ection is an effect of deflection of high- 
energy charged particles upon their over-harrier (non- 
channeled) passage through a planarly oriented bent 
crystal. The effect arises when the crystal bending radius 
R greatly exceeds the critical value Rc- That condition 
is the same as the Tsyganovs one for the possibility of 
channeling in a bent crystal py], but the particle motion 
regime in the crystal yet depends on the particle entry 
angle relative to the active atomic planes. When this an- 
gle is much larger than the critical value 9c, then, moving 
in the continuous potential of bent planes, conserving the 
particle full transverse (radial) energy, the particles are 
rarely captured into channels (via incoherent scattering 
on atomic electrons and nuclei), and are mostly deflected 
elastically through the volume reflection mechanism. Cu- 
riously, the latter deflection proceeds to the side opposite 
to that of the crystal bending; the value of the deflection 
angle is of the order of critical angle Oc- Furthermore, 
the particle beam after deflection remains fairly well col- 
limated, i.e., its angular dispersion keeps much smaller 
than the mean deflection angle. That phenomenon was 
discovered in numerical simulations two decades ago [2] 
and recently verified experimentally [1, Q • Nowadays it 
is considered to be an option for beam coUimation and 
partial extraction at ultrarelativistic charged-particle ac- 
celerators & 



To a good accuracy, the particle dynamics in the vol- 
ume reflection problem is classical @ and reduces to 
classical particle motion in a cylindrically symmetrical 
continuous potential of bent atomic planes. Therewith, 
granted the angular momentum conservation relative to 
the active crystallographic plane bending axis, the final 
deflection angle is expressible in the standard way as an 



'Electronic address: ' bon@kipt.kharkov.ua| 



integral over the radial coordinate from an inverse square 
root function involving the potential [see Eq. (|55bl) be- 
low]. That approximation served as a starting point for 
a number of numerical studies 0, Q . 

Although the described computational problem seems 
to be sufficiently simple, it is aggravated by the presence 
of several parameters: the ratio R/ Rc, initial and final 
particle variables (impact parameter, the angles of in- 
cidence and deflection). The dependencies on all those 
parameters involve singularities, which in general are bet- 
ter dealt with by analytic techniques than by numerical 
ones. Besides that, as long as for practice mostly inter- 
esting is the case i? ^ i?c, it would be instructive to eval- 
uate the asymptotic behavior of all relevant observables 
in the formal limit R/ Rc — ^ oo, including next-to-leading 
order corrections in the small parameter Rc/ R. But since 
volume reflection depends on the particle dynamics not 
in one but in several inter-planar intervals, for feasibility 
of its global analytic description one rather needs some 
model. 

A valuable opportunity for realistic model building is 
that the inter-planar potential in a silicon crystal, at least 
in the orientation (110), is fairly close to parabolic s hap e 
over the entire inter-planar interval (see, e. g., jlC|). 
A parabolic (harmonic) potential, i.e., a linear oscilla- 
tor, permits a simple solution for the particle trajectory 
within a single interplanar interval. The next problem 
is to connect solutions on the boundaries of the adjacent 
intervals. It may appear nontrivial, but it is feasible to 
do that transitively, i.e., simultaneously for an arbitrary 
number of the adjacent intervals. Thereby we obtain a 
completely solvable model capturing basic features of the 
volume re?ection, except the effects of incoherent multi- 
ple scattering. Moreover, we are able to derive not only 
the de?ection angle, but also an expression for the whole 
trajectory, which further on may be used for description 
of inelastic processes, such as volume capture or electro- 
magnetic radiation. 

In the present work we will deliver a solution for the 
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posed model problem. The plan of the article is as 
follows. In Sec. |lT] we describe the procedure of solu- 
tion connection between adjacent inter-planar intervals, 
demonstrating that the problem reduces to elementary 
trigonometry. The particle trajectory is expressed as an 
explicit function of inter-planar interval order number 
(not involving a recursive procedure) , for arbitrary ratio 
R/Rc- In Sec. IIIIl from the obtained solution for trajec- 
tory, we derive the particle final deflection angle, which 
comes as a sum of inverse trigonometric (for positively 
charged particles) or hyperbolic (for negatively charged 
particles) functions. In Sec. IIVI we scrutinize the limit 
i? 3> i?c, most interesting in relation to volume reflec- 
tion and practical applications, first for positively, then 
for negatively charged particles. In the generic expres- 
sion for the ?nal deflection angle, we flnd a possibility 
to replace the sums involved by integrals (via the Euler- 
Maclaurin formula), and do the latter ones in a closed 
form. As a result, we arrive at sufficiently simple asymp- 
totic formulas for the de?ection angle dependence on all 
the variables. The impact parameters are thereupon an- 
alytically averaged over, and the experimentally observ- 
able scattering differential cross section is obtained for 
positive and negative particles. In Sec. |V] we examine 
the opposite limit i? <C -Rc- In Sec. |Vl] we provide es- 
timates of optimal crystal and initial beam parameters 
for beam complete dcfiection or for experimental investi- 
gation of the final beam profile features. A summary is 
given in Sec. IVIII 



II. PARTICLE TRAJECTORY IN A BENT 
CRYSTAL 

A. Initial conditions 

The usual geometry of experiments on volume refiec- 
tion iniplies sending a charged particle beam normally to 
a thin [20|, weakly bent crystal plate. The practically un- 
avoidable slight curvature of the crystal boundary thereat 
is of minor consequence, since the main contribution to 
the particle refiection angle comes from a vicinity of some 
point in the depth of the crystal. For definiteness and to 
establish an easy connection with the particle impact pa- 
rameter in the initial (perfectly parallel [2l|) beam, let us 
consider a particle incident along the z axis on a crystal 
whose front face is a perfect plane, located at z = 0. As 
for the crystal rear face, for our purposes in this paper 
we may leave it unspecified at all, as if the crystal was in- 
finitely thick, but transparent. Then, let (0 < ^ 1) 
be the angle of inclination of crystalline planes to the z- 
axis at the crystal front face (see Fig. [Ij, and let a;-axis, 
perpendicular to Oz, point in the direction of the crys- 
tal bend. Moving at small angles to the crystal planes, 
the particle interacts most strongly (coherently) with the 
averaged, so-called continuous inter-planar potential Q, 
which induces a force with dominant a;-component (yet 
slowly dependent on z) [22l|. Along the y coordinate there 
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FIG. 1: Coordinates describing the bent crystal geometry (cir- 
cle segments are the bent atomic planes) and the particle en- 
trance to the crystal (thick arrow). Not to scale. For details 
see text. 



is a translational invariance, ensuring the particle mo- 
mentum y-component conservation. 

In monocrystals of not too heavy chemical elements, in 
particular for silicon (lattice of diamond- type) , oriented 
by its (110) plane close to the direction of the beam, 
the continuous potential in each inter-planar interval may 
closely be approximated by a quadratic function, with an 
accuracy < 20% [1^. That entails a linear equation of 
motion for the classical [13] ultra-relativistic [l^ particle: 
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(c = 1, small angle motion). 



(1) 



(2) 



where is the midpoint of the inter-planar interval, d 
the inter-planar distance, E the particle energy, and i^max 
the force acting on the particle at the edge of the inter- 
planar interval x — xq — — |. For positively charged 
particles i^max is positive, whereas for negatively charged 
particles it is negative. Note that the force and the par- 
ticle energy enter equation ^ only through the ratio 
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known as the Tsyganov critical radius p|. The natural 
time unit in channeling-related phenomena is 
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(27rr has the meaning of positively charged particle chan- 
neling period, although herein we deal not with channel- 
ing but with an over-barrier motion) . 

One and the only consequence of the crystal bending 
is that xq in Eq. ([1]) acquires dependence on the longi- 
tudinal coordinate z, which for ultra-relativistic motion 
under small angles to Oz may be equated to the current 
time t: 



xq = xq{z k, t) (the crystal bend function) 



(5) 
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In application to volume reflection, we are interested in 
the uniform bending of the crystal, at which Xo(t) de- 
scribes a circular arc of a small opening angle. That small 
arc may equally well be approximated by a parabola, and 
hence xo{t) is determined by the equation 

xo{t) = —Oat -\ (uniformly bent crystal) , (6) 

2R 

where R is the atomic plane bending radius (without the 
loss of generality one may let a;o(0) = - see Fig. [T]). 
Inserting (jH]) to ([T]), and implementing Q, we get the 
particle equation of motion in the first inter-planar inter- 
val: 
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pos. charged particles 1 
neg. charged particles J ' 
(7) 



Initial conditions for x{t) stand as 

x{0) = h, 
x{Q) = 0, 

where h, restricted to the interval 



(8) 
(9) 



(Here and henceforth upper signs and figures refer to pos- 
itively charged particles, and lower ones - to negatively 
charged particles.) Matching initial conditions (|14m5|) 
allows one to determine constants and ip^'- 
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arcsm 
arsinh | 



(17) 



(18) 



(the sign of Ao must be chosen positive in order that f(0) 
at ^0 > 0, according to (ITSI) . was negative). 

Further on, solution ([TS]) is to be connected with so- 
lutions in the subsequent inter-planar intervals. Impor- 
tantly, since these solutions have to be connected at a 
definite r, the condition of connection will not depend on 
the current phase of the harmonic motion, such as lpq in 
Eq. (|16p . as we are going to show. 



B. Connection of solutions through interval 
borders 



d , d 
- < 6 < -, 
2 - - 2 



(10) 



is the impact parameter measured from the middle of the 
interval 

The equations of motion further simplify in terms of 
the "subtracted radius" variable 



r{t) = -x{t) - 6'o< + 



becoming 



where 



5^r 



2R' 



d d 
- < r < - 
2 - - 2 



(11) 
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Moving along trajectory p^ . the particle will cross the 



next inter-planar interval border r 



I'l I arcsm | f ^ 



at an instant 
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^0 



^0 



(19) 



(inferred from (IT51) by letting r = — | and solving for t). 
At this instant, the equation of particle motion turns to 



dTid + r) 



M d 

m < r < 

2 - - 2 



(20) 



That is the same harmonic oscillator, only with an altered 
equilibrium position, and general solution of (j20p may be 
written as 



(13) 



ri{t) ^±S-d-Ai 



sm 
sinh 



t 



po + ^fi]- (21 



Thus, ±5 is the spatial shift of the oscillator equilibrium 
position due to the crystal bend, i. e., due to the cen- 
trifugal force, which in the present small-angle approxi- 
mation, presuming condition r <C i? within the weakly 
bent crystal, is treated as virtually independent of the 
subtracted radius r (cf. 0). For r(i), the initial condi- 
tions ([HE]) translate to 



r(0) = -6, 
Generic solution of Eq. ([T2|) reads 



- + fo 

T 



(14) 
(15) 



(16) 



Values of new constants Ai, Aipi are now to be deter- 
mined from the continuity of r{t) and f(<) at the border 
point r — That can be done without formally solving 
the system of two equations. First, compare two integrals 
of motion 



Af = T^rl ± {±S - d - nY 
Al = T^rl±{±S~rof 



(22) 
(23) 



(related to transverse energy) in their common point, 
where ro — ri = — |, fo = f'l- Subtracting (|23l) from 
I, one gets A^ = A"^ — 26d, i. e. 



Al 



Al - 25d. 



(24) 



FIG. 2: Solid curve - the phase-space (the subtracted radius 
r vs. the radial velocity f) trajectory for positively charged 
particles, under condition 5 < |. Dashed vertical lines signify 
the positions of bent atomic planes (definite r). Thick dots 
indicate centers of the trajectory circular segments. Vertex 
angle Ayi (and similarly all other Aipn) may be interpreted 
as a geometric sum of vertex angles in a pair of right triangles 
having a common cathetus, and with the second catheti equal 
I + (5, f — 5, and the hypotenuses Aq, Ai. 



Thereupon, the phase shift A</?i is sought from a condi- 
tion ri(ti) = — |. One finds: 

/\ — { ^-I'csin 1 f i f arcsin 1 f T , . 
^ I arsinh J Aq \ arsinh ] A\ 

(for a geometric interpretation of this relation for positive 
particles - see Fig. [2]). As we had expected, neither A\, 
nor A(/3i depends on Lp^. 

At each next border connection of the solutions is im- 
plemented in exactly the same way. Writing in the rt-th 
interval 

TniPi = ±(5 - - A„| I ^~ + + X! ^"^"i j ' 

(26) 

( d ^ d ^ \ 

[~--nd<r<--nd, t,, < t < t^+i I 

the generic amplitude is found as 



A„ ^J Al_^ - 25d ^ ^ Al - 2n5d 

= ^jT^e^±ib±Sy -2nSd. (27) 

and the phase shift is deduced to be 

. f arcsin 1 | =•= (5 ( arcsin 1 | T <^ , . 

^'^"--(arsinh/A:::;- I arsinh/^' ^^8) 

where the amplitudes in the denominators must be 
treated as already known, by (|27p . The instants of border 



FIG. 3: The same as Fig. [21 but for negatively charged parti- 
cles, assuming condition 5 < f (see Sec. IIIIBj) . 



passage can also be evaluated: 

tn _ J arcsin 1 f i J arcsin 1 i T <^ 

— - 2^ [ arsinh / A^ ~^ arsinh j A„ 

One caution is that amplitudes An should not be re- 
garded as a measure of the particle spatial wiggling in 
each interval. As Figs. [21 [3] indicate, the trajectory swing- 
ing enhances as the particle penetrates deeper into the 
crystal, whereas amplitudes An, to the contrary, decrease. 
There is no contradiction hereat because for most of the 
intervals traversed the intra-channel oscillation period 
2ttt is much greater than the time of particle passage 
across the interval, thus the particle is far from mak- 
ing a full oscillation in each interval, anyway. In fact, 
the lower is the amplitude An compared to the interval 
length (along the particle motion direction), the stronger 
warp of the trajectory on this interval may occur (see 
Figs, [a [3]). 



III. PARTICLE REFLECTION 
A. Reflection conditions for positive particles 

It is clear that the decrease of amplitudes (|27p can 
not continue indefinitely, because eventually arguments 
of the arcsines in ([28]) shall exceed unity (that happens 
sooner than the radicand in ([?7)) becomes negative). This 
merely signals that the particle can not reach the next 
inter-planar interval. The particle will continue the har- 
monic motion until it hits the previous interval, then pro- 
ceeds moving outwards in the radial variable in the same 
but reverse way, and on the exit from the crystal it will 
emerge as a deflected beam. 

Let us evaluate the order number Umax of the reflec- 
tion interval. If for some n inequality ^ + 6 < An-i is 
met, then it also entails ^ — S < An, so arguments of all 
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the arcsines in ((29|) are less than unity. So, is the 

largest integer yet allowing for ^ — S < A (+) . Through 

(l27)l , that condition determines the reflection interval or- 
der number: 



,(+) 



.2/12 



2Sd 



(30) 



where the lower-corner brackets [. . .J designate the inte- 
ger part of a number ( [aj < a). If t9o ^ d, variation of 
Timaxib) is much smaller than its mean value 



2Sd 



2d 



(31) 



see 



(essentially valid for negative particles as well 
Eq. ([?T|) below). It is instructive to notice that 



- - (J < A (+) < -+S. 

2 - "max 2 



Toward the volume reflection problem, we are inter- 
ested in finding the total reflection angle ^rofl, half of 
which, by symmetry reasons, amounts to deflection an- 
gle of the trajectory in the reflection point t — t^cfi in 
which 



r(irofi) 



0, 



(32) 



< 1), one recasts (p4)) as 



ticR TT 

— — ~ arcsm 

r 2 



b + S 



^(|-<5)' + 2(t.(+)+4il) 5d 



E 



V ^{l+5)' + 2{.i+)+n) 



Sd 



- arcsm - 



^(f -<5)'+2(z/(+)+n) Sd^ 



(36a) 



Equivalent ly, using the identity arcsin -^^^ — arccoty^, 
one can write 

b + S 



trcfi TT 

= arcsm 

T 2 



^(1 -<5)' + 2(j.(+)+nLU) Sd 

^ / j2{i^i+)+n)Sd J2{iy(+)+n)Sd\ 
+ > arccot ^ ^ , l-arccot ^ ^ , — . 

(36b) 

The physical meaning of parameter v^'^'^ is clear from 
Fig. [2] It represents the kinetic transverse energy at the 
last atomic plane before the reflection, divided by the 
centrifugal potential difference between the neighboring 
atomic planes. 



i^rofi ^ i(ircfl) = -00 + S^^ (33) 

(more exactly - see Eq. ((52)) below) [2^. To evaluate 
the right-hand side of ([33]) . one only needs to know the 
value of irofl. The latter is found from solving equation 
r(ircfl) = 0: 



r 



r 2 

; . , arcsm - , 

2 ^ \ An-l 



arcsm ■ 



b + S 

AT' 



til 

An 



(34) 



The largest contribution to the emerging sum comes from 
the terms n ~ Umax (where denominators An are small- 
est), so it may be more convenient to revert here the 
summation order. Introducing a useful parameter 
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9i + {b+sr-{i-sy 

2Sd 



(35) 



with braces {. . .}f to indicate the fractional part (0 < 



B. Negative particles 

In contrast to trigonometric arc-sine, hyperbolic arc- 
sine function exists at any value of its argument. There- 
fore, expression p6)) for negatively charged particle tra- 
jectories holds until the radicand in the motion amplitude 
An given by (P7|) becomes negative. The first interval at 
which that happens will be called "infiection" one. Its 
order number is inferred to be 



ninfl 



.2/12 



[b-sy 



2Sd 



+ 1. 



(37) 



In the inflection interval the amplitude An^^^^ calculated 
by the formula (j27p would be imaginary. That implies 
that the r{t) dependence now is to be described by a 
hyperbolic cosine rather than a sine (hence the term "in- 
flection" ) . Matching the amplitude and the phase of the 
hyperbolic cosine with solution (|26p for the preceding 
n ^ riina - 1 gives 

A «i„fi-i 

''ni„fl(^) = -<5-"inflrf+P^ni„fil cosh -+(y3o+ ^ Aipm 

\ m=l 

^ - S - +s\ 

— arsinh— 2 arcoshv^ 
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"infl- 

d , d , 



I^Mnfl I / 



(38) 
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where 



with 



2Sd 



(39) 



(40) 



Therethrough, using Eq. ((33)) . results the deflection an- 
gle. 



Otherwise, i. e. if ^ holds (e. g., if f > 3 > /), 
in all the subsequent intervals after (|57)) the trajectory 
must also express through hyperbolic cosines: 



Since ^ + 6 > V26d, the argument of arcosh in ([38|) is 
always > 1. 

Next, a question arises, at which condition the tra- 
jectory ([55)1 can actually reach the next interval, i. e. 



''"infi (t) can descend to value ninfld— |. Since by Eq. ([iU]). 
r„jjjj(t) > —(5— riinnd + |A„j^jj|, that would require 



S~->\A„ 



(41) 



Substituting here (P^)) . and solving with respect to the 
ratio 4, one may present (j4ip in form 



2R 



(42) 



with 



/(^(-^ - ^ - ^'"^ + v'(2-'^(-')(l-^(-0- (43) 

Function /(i/*^^)) decreases monotonously (almost lin- 
early) from /(O) = I a/2 « 2.9 to /(I) = i. 

In the simplest case illustrated in Fig. [31 when condi- 
tion dH]) is violated (e. g., if ^ < i < /), ([271) must be 
the last interval reached by the particle, its order number 
being 
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25d 



1 



(44) 



Expressing ticfi from equation r(trcfi) = then gives 



^rcfl , f + 

= arcosh-^; 

r L4 



arsinh- 
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arsinh - 



(45) 



Reversal of the summation order here leads to an expres- 
sion 
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■ arsinh — ^ 
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(46) 
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r„(t) = — (5 — nd + |A„| cosh j — h (^so + ^ A 



m— 1 



- - 5 - -I- 15 

-arsinh— 2 arcosh-^ 



^ninfi — 1 

+ 5 



I^Minfl I 



Earcosh-T- — r -f arcosh , , 

\Am\ ^ \A, 

- - - nd < r„ < - - nd, n> ninfl 



4 + '^' 



(47) 



with amplitudes An still given by Eq. (l27t . Sequence 
()47p may continue as long as the arguments of all arcosh 
exceed unity, i. e. as long as 



d 



(48) 



(which is equivalent to — | -f ^ > |A„_i|). Inserting here 
(P7)l . one ultimately infers the value of the reflection in- 
terval order number: 



,(-) ^ 



.2/12 



2 _ (5-^)2+^1+^ 
2M 



1 > fi^'-' 



,(49) 

The above expression is similar to Eq. pop for positively 
charged particles. As one might expect, in the high- 
energy limit 5 ^ d,b values ni^^x and rtinax coincide and 
do not depend on the particle energy. 



Expressing tj-cfl from r(ircfi) = and Eq. (|47p in this 
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case gives 



C. Thick crystal limit (isolation of volume effects) 



tied . , b - 5 

= — arsmn — : — 

T Ao 
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- arsinh^— h arsinh^— — . 
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arcosh- 



— — arsinh- 



\A,n\ 
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(50a) 
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- arsinh — ^ , 

n=o ^2Sd (iy(-) + n) 

"L";'x-nintl + l + (5 

arcosh ^ 



(-)) 



arcosh 



y^2M (n-U-)) 



(50b) 



Actually, equations (|50bp can be used not only under 
condition (j42|) . but also at any ratio ^, provided that in 

capacity of n,nax one uses expression 



,(-) 



ib-sr + ii-sfQ{i-f{.i-^)) 

2Sd 



(with 0{v) the Heavyside unit-step function) unifying 
(|44|) and (j49| . The universally valid formula (jSTj) may be 
convenient when iicfi is evaluated with the aid of com- 
puter for widely changing values of particle energy or 
crystal bending radius. 

As for the physical meaning of at ^ < /, i. e. 

when inflection interval is also that of reflection, Fig. [3] 
illustrates that the meaning of is similar to that of 
It is the (appropriately rescaled) kinetic transverse 
energy upon the particle entrance to the reflection inter- 
planar interval, only the interval boundary now is not the 
atomic plane but the last potential maximum passed. In 
case if I > /, i/'^^ does not characterize the reflection in- 
terval, and vice versa, the kinetic energy in the reflection 
interval is not closely related with i^'^^. 

The obtained expressions (|26ll29p for the trajectory 
and (|36l I50bp for its reflection point allow evaluating all 
the observables relevant to the particle passage. In the 
present paper, we will be interested only in the final angle 
of elastic reflection. 



Formulas (|36[ I50b| . in principle, contain dependencies 
both on volume and on boundary effects. In most practi- 
cal cases, the deflecting crystal may be regarded as thick, 
whence boundary effects are expected to turn negligible. 
An increase of the crystal thickness, or more precisely of 
the distance between the crystal boundary and the vol- 
ume reflection point, may be thought of as an increase 
of the particle incidence angle 6*0 (see Fig. [1]). Then, it 
suffices to consider the limit 



2 hm 

eo/Bc^< 



-^rcfl(^o) 



(52) 



With function (1551) . or (|50b|) . such a limit must always 
be finite: indeed, at large ^max the sum over n grows like 
the corresponding integral, whose asymptotic behavior 
straightforwardly evaluates as 



T 



dn 



1 + 6 



2 " 



^J2{v^+) + n)5d v/2(ty(+) +n)5d 



2n„ 



This leading asymptotic behavior cancels exactly the first 
term in (1521) . whilst calculation of the finite remainder 
requires a more accurate evaluation of the sum, which 
will be our task in the next section (in application to the 
limit i? > i?c). 

In general, it must be noted that function 
^'v.r.(''', <5, d, 6), being a dimensionless function of 4 
dimensional variables, may depend only on their 3 
dimensionless ratios - say, d/r, 6/d, b/d. At that, the 
last ratio is always '-^ 1. The first one amounts to 




261, 



(53) 



where 9c is the Lindhard critical angle [8|; so, it is always 
small, once we are in a high-energy regime. As for the 
ratio 



25_ 
~d 



Rc 
li' 



(54) 



it may be either large or small depending on the particle 
energy and the crystal bending radius. The regime of 
particle passage through the crystal is determined solely 
by ratio dM]) . 

To gain a general impression of functional dependen- 
cies involved, and to test our generic formulas ([551 ISObp . 
let us view the dependence 0v.r.(&) for different values of 
R/Rc- Fig. |4] shows this dependencies for positive and 
for negative particles. They are in fair agreement with 
Figs. 6 and 8 of [ll|. But we will pay more attention to 
interpretation of the features observed in the figures: 

1. The origin of the recurrent structure in variable b 
with a tapering period is, obviously, due to O^.r. de- 
pendence on b through z/^^' alone (see Eqs. p5l) . 
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2Rld 



potential V{r) j28| : 





FIG. 4: The volume reflection angle Sv.r. as a function of 
impact parameter h, for R/Rc = 1/3,2,10,40. Top panel - 
for positively charged particles, bottom panel - for negatively 
charged particles. Beyond the shown unit interval the picture 
repeats periodically. 



(PO))). insofar as involving an operation of 

fractional part, is a periodic function of 
which in the interval — ^ < 6 < 4 makes 



2Sd 



4<b<i 

2^ ^ 1 periods. Also, since t^'^*-' is an even func- 
tion of 6 ± (5, the particle deflection angle is a sym- 
metric function of b with respect to point b = —S 
for positively charged particles, and with respect to 
6 = i5 for negative particles. 



Another feature of 0v.i.{b) dependencie(s) is that, 
for negatively charged particles, the reflection an- 
gle blows up (formally) to -f oo at certain values of 
impact parameters. Physically, that corresponds to 
close matching of the particle transverse energy to 
the height of a (locally parabolic) effective poten- 
tial barrier - the situation known as orbiting (see 
[l3|) [131 • The asymptotics of the divergences is 
logarithmic fl^, as follows from the general inte- 
gral expression of the deflection angle in a central 



2M 



1 



dr/{R + rf 



^{E-V{r)f-NP/{R+rf 
dr 



rbit. 



4^ 



si _L V.ii{b)-V^ii{r) 

dr 



10. , 

4 ' 2E 



I IP I ' I,?' '"saddle/ 



TT (55a) 
(55b) 
- . . . (55c) 



4(5 J In^ -t- ... (A£;j_ cx ^ +0) 

5 1^1^^ + - (A-Bi cx - 1 ^ -0). 

(55d) 



Here 



M~{R-b)[E- V{b)] (1 - 61^/2) 

is the particle angular momentum relative to the 
crystal bend axis. 



V,s{r) = V{r) - E 



R 



(56) 



is the effective potential including the centrifugal 
energy, rgaddic - the position of maximum of the 
effective potential barrier whose height in the case 
of orbiting happens to be close to the particle en- 
ergy, and Ai?_L - the transverse energy variation 
relative to the height of the effective potential bar- 
rier. The factor i in the lSE± < alternative of 
Eq. (|55d[) arises because the integration in (|55cl) is 
then carried out only over the one-sided neighbor- 
hood of Tsaddic where the radicand stays positive 
(see Fig. El. 

3. It must be noticed that for negative particles func- 
tion 0V.I.Q}) li^s smooth minima, which must corre- 
spond to caustics, i. e., to rainbow scattering [12]. 

4. In contrast, for positive particles the potential in 
its maximum is not differentiable, excluding both 
orbiting and rainbow scattering for this case. With 
some smearing of the potential around the atomic 
planes, these affects, of course, re-appear. 

The proper question is whether it is possible to de- 
rive at least the particle final deflection angle (related 
with ^^) from the more conventional integral represen- 
tation approach 0, Q . In that approach momentum and 
transverse energy conservation laws are incorporated au- 
tomatically, so there is no need to connect trajectories on 
the interval borders. Indeed, specializing in (I56[ I55bl) 



V{r) = F^^^d 

(at r < n > 0) 









r 


1 


71^ , 


n = 




rd' 


" 2 



(57) 
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eff 




'"saddle 



to determine the numerical value of the limiting ratio 
limfl./fl.^_j.oo Secondly, it is desirable to determine 

the final beam shape and quantify its angular width as a 
function of R/Rc- That will be our aim for the present 
section. The treatment is somewhat different for the 
cases of positively and negatively charged particles, be- 
cause of the difference between the functional form of 
initial Eqs. §^ and (l50b)) . 



FIG. 5: The relation between the particle transverse energy 
and the effective potential energy T4ff (including the centrifu- 
gal potential) under the conditions of negatively charged par- 
ticle orbiting in a bent crystal. The leading logarithmic con- 
tribution to integral H55b|) comes from the vicinity of point 
''saddle - the coordinate of the effective potential maximum to 
which the particle transverse energy happens to be close. In 
the case AEj_ > the particle sweeps the two-sided neigh- 
borhood of Tsaddic, whereas at A_Bj_ < - its one-sided neigh- 
borhood only. 



and using basic integral 

d/2-nd 



dr 



-d/2-nd ^T^el + 62 + 2(5(6 + r) - (r + ndy 



arcsm ■ 



-|- arcsin 



y^T^e'^ + (6 + (5)2 - 2n6d 

l + s 

^Jt'^BI + (6 + 5)2 - 2n6d 



(58) 



for positively charged particles, and a similar one for neg- 
ative particles, we reproduce the inverse trigonometric 
and hyperbolic functions encountered in ([Ml H51 I50al) . 
But the integral representation approach wouldn't give 
us explicit trajectories r(i) (rather, t(r), to be solved for 
r), and the geometric interpretation (Figs. [2l [3]). 

On the other hand, from the integral representation 
for the final angle we might derive the result in form 
of a sum of analytic functions also for a more compli- 
cated parametrization of the inter-planar potential - e. g., 
adding thereto a term proportional to r^. Then instead 
of arcsines one would get elliptic functions. But it is the 
simplicity of functions in the sum that permits us, in the 
important limit R ^ Rc, when the number of terms in 
the sum gets large, to replace the sums by integrals and 
do them in closed form. In this sense, analytic investiga- 
tion only begins here. 



IV. VOLUME REFLECTION REGIME 
(MODERATELY HIGH ENERGIES, R » Rc) 

As we had mentioned in the Introduction, and as Fig.|4] 
does confirm, under the condition R ^ Rc, i. e. 26 ^ d, 
the particle deflection angles depends weakly on the im- 
pact parameter. So, it is interesting, in the first place, 



A. Positive particles 



For positively charged particles, in the considered lim- 
iting case 6 d, say, quantity (| ± (5) +2 -t- n) Sd 
entering the denominators in (j36a|) varies relatively little 
as n changes from n to n + 1. Thus it must be possi- 
ble to replace the summation in (j36p by integration. The 
proper mathematical tool for that is the Euler-Maclaurin 
formula (see, e. g., [l3|) which reads 



E/H=2^(0)+/ dnf{n) + - 



f{N) + 



dn 



(59) 

Employing this formula for approximation of each of the 
sums in (j36bp (which is somewhat more convenient than 
Eq. ((Ma|) ). one gets [H 



arccot 



^25d {vi+) + n) 



n=0 



i±5 



1 y/25dvW 
— -arccot — -1 

2 !k±5 



,(+) 



dn arccot 



\/2Mn 



^ ^25d{n^^l + u^+)-l^ 



-arccot 
2 



±5 




where we had estimated, for all n, 



dn 



arccot 



±5 




(60) 



The two end-point contributions in the third line of (|60p 
are small as O {S/d) relative to the integral, but still they 
need to be kept if we wish to describe not only the mean 
deflection, but also the scattered beam shape. 

Taking the indefinite integral in Eq. (pO)) by parts 
J dn arccot Van = — [ (l-f an) arccot \/an -I- Van] , 
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one brings (|60| to the form 



arccot 



^25d{vW + n) 



n=0 



±5 



25d 



±sy 



,(+) + - 



2 I max 



•arccot 



^2(5^ (nLU + - l) 



l±5 



l±6 



arccot - 



±(5 



±6 



1 V2^ 1 y2^<4iU.W-l) 



--arccot 
2 



, r -arccot - 

^±5 2 




(61) 



In the hmit n^^^ — oo, with the use of asymptotic ex- 
pansion arccoty^ ^ f ^ + ^ {v^^^) : expression ([6T] 
reduces to 



^26d + n) 






I I I I I 



FIG. 6: Solid curve - positively charged particle reflection 
angle Sv.r. versus impact parameter b, for R/Rc = 25. Dashed 
curve - approximation (|63|) . The central segment of the curves 
is strongly So-dependent. 



ultimately arrive at result 



..r.(6) 




l-l|(l-.W(6) 



R. 



3/2' 



i?3/2 



(63) 



(with iy(+)(fe) being given by Eq. ([35])). 

Comparison of approximation (|63p with the exact re- 
sult ([36t is shown in Fig. |6l (Actually, the given ap- 
proximation appears to be numerically accurate starting 
from R/Rc ^ 5). From the figure (or Eq. one con- 

cludes that in the first approximation all the particles 
are deflected to the same angle ss — f ^c- There is also 
some dispersal of the scattering angles, depending on the 
particle impact parameter, of the full width 



27TRr 



2ttS 



R 



(posit, charged particles) 
. 

The observable quantity, however, is not the indicatrix 
^v.r.(^) but the scattering differential cross-section (fi- 
nal particle flux averaged over the impact parameters b) 
as a function of the scattering angle ^v.r.- Therefore, it 
is desirable to reconstruct the latter dependence issuing 
from the first. That does not pose any principal problem, 
granted the linearity of dependence 6'v.r. (i^^^-*)- 



V2SdM+T , V2Sdu(+1 



the 



Here, one notices that terms — , -i ttt- 

brackets in ((62|) cancel. Further on, inserting (|62|) to 
(|36b[) and this to ([5^ . we witness the anticipated can- 
celation of the large terms — + -^^-^"-max — 0, and 



1. Differential cross-section 

Turning to evaluation of the differential cross-section, 
one encounters a certain complication: the 6-dependent 
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quantity u^~^'> in (|63p also contains dependence on ft). 
The fact of residual 6*0 -dependence was noticed in pl| . 
To some degree, it conflicts with our initial assumption 
about the boundary condition vanishing influence in the 
limit of large 9o/0c- We can not revoke it at the present 
stage, since in equation ()52p we had already incorpo- 
rated the facilitating assumption of the trajectory sym- 
metry with respect to point ii-cfl. Obviously, the sen- 
sitivity to the boundary conditions in general destroys 
such a symmetry. Furthermore, ^refl might as well con- 
tain a dependence on the particle exit angle relative to 
the atomic planes, which we did not even take trouble to 
specify. Altogether, that may rise a suspicion that the 
^-dependent correction obtained in (|63|) is unreliable for 
evaluating the differential cross-section. Fortunately, the 
impediment is not fatal and can be overcome within the 
present framework. In principle, the differential cross- 
section sensitivity to 9q is attenuated with the increase 
of R/Rc, but more importantly, we will prove that upon 
averaging over a tiny interval of 9q this dependence is 
eliminated completely. 

To begin with, the differential cross-section involves 
only a derivative of function 0v.r.(&): 



dX 



E 



1 



M6'v.r.M|f,^f,„(e_) 



di^W/db\ 



(65) 



b=b„(iy( + )(6lv.r.)) 



where bm{d) is the set of all the roots of equation 
dv.r.{b) = 6 belonging to the interval -| < 6 < |. Now, 
at R/Rc^ 1 the number of roots bm of equation ^v.r. = S 
is large, and so, in general, they are densely distributed 
over the finite definition interval — | < 6 < | . It appears 
that the root distribution density is just proportional to 
the derivative in the denominator of (1651) (the formal 
demonstration is relegated to the Appendix). Therefore, 
the sum appearing in (j65|) approximately equals to just 
the b variation interval length, i. e., d. However, the 
relation expected thereby. 



dX 
ddv.T 



Rd 



2-kOcRc 



(66) 



does not yet hold uniformly in 6, and hence in ^v.r.- For 
instance, in a neighborhood of point b = —5 we have 
in the denominator of (|65|) dv^^'> /db (see Fig. [6]), 
so there the differential cross-section blows up above the 
plateau ([TO)) (see Appendix). But the latter peak position 
on the ^v.r. axis depends sharply on the value of 9q and 
hence is essentially "random" , needing to be averaged 
over. 

Indeed, one notices that the dependence of 1^*^+' on 
^0 is quadratic, so a situation is possible when the inci- 
dent particle beam divergence is smaller than the angular 
spread acquired in the crystal: 



9, dA 
d dft,. 



R 



(1- 



FIG. 7: Asymptotic (at R ^ Re) behavior of the So-averaged 
differential cross-section for positively charged particle scat- 
tering (Eq. (|71|) V The area under the rectangular curve is 
unity, representing the total probability. In higher orders in 
Rc/R the distribution edges must smear out (see discussion 
in the text). 



but at the same time, the indeterminance of ~ ^max 



is greater than unity: 



A 



2Sd 



Sd 



A6'o > 1. 



(68) 



Together, Eqs. (|571 [55]) may be viewed as a double in- 
equality 

— ^ < A0O < — (00 -averaging). (69) 
Here the sufficient gap exists provided 



(70) 



This is basically the same condition that we had assumed 
in writing Eq. ([5^ . thus for derivation of a 6'0-averaged 
differential cross-section we can safely rely on Eq. (|63|) . 

Ultimately, we can make a specific statement that un- 
der conditions (1691) . upon 0o-averaging, the differential 
cross-section equals to constant (|66p over an interval 
where roots bm exist. This interval is determined in the 
Appendix (Eq. (theta-interval)). So, the 6'0-averaged dif- 
ferential cross-section (the final beam profile) is described 
by a simple rectangular function 



d9^ 



Rd 



^TrOcRc 
n 



TT 

2 



AR, 
~R 



(71) 



(see Fig. III). 



2. Comparison with experiment 



(67) 



The deflection angle mean value is least affected by 
multiple scattering, and thus, may be directly compared 
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with the experiment. From (|63p we obviously infer 



= -^eA 1- 



2R, 



d 
2R, 



1 - 



2R^ 



(72) 

The property of ()72p is the Unearity of dependence on 
the crystal curvature R~^] the linear kind of depen- 
dence was indeed noticed in CERN experiments with 
E = 400GeVGeV/cm 

To make quantitative comparison with the experiment, 
one needs to specify the potential strength in our model. 
In reality, the Si (110) inter-planar potential is character- 
ized by 2 parameters: F^ax ~ 6GeV/cm (usually used 
for evaluation of for channeling processes) and the 
well depth Vq — 22.7 eV (usually used for evaluation of 
the critical angle 6c for volume reflection). The relation 
Vb = -j^lnaxrf implied by a quadratic potential model 
only holds with accuracy ~ 20%: 

22.7 eV = 7^ J^maxd - 28.8 eV 
(d= 1.92A for Si (110)). 

If we evaluate Rc in (17^ as Rc = E/F^a.^ ~ 0.67 m, it 
will produce too large |(0v.r.)|- But evaluating both Rc 
and 6'c as i?c = ^ « 0.85 m, Oc = ^JWq/E « 11 ^rad, 
and substituting to Eq. (|72t , we get a satisfactory agree- 
ment with the experiment (see Fig.[8|). 

As for the obtained rectangular profile shape, it is more 
sensitive to multiple scattering, and was not yet probed 
by experiments (the optimal experimental conditions will 
be specified in Sec. IVip . But we can compare our profile 
with the available numerical simulation results using a 
realistic, smeared potential, without multiple scattering: 
Fig. 6. In that case, the positive particle profile 
shows indeed a signature of flattening ( "shoulder" ) but 
near its edges the distribution behaves differently, ex- 
hibiting a subtle divergence (rainbow) at the outer edge, 
and decreasing continuously on the inward side. So, 
for positively charged particles our simplified model of 
parabolic inter-planar potential describes the final beam 
profile only quantitatively, though it is able to predict 
the distribution width and mean value. 

Next, turning to the negative particle refiection prob- 
lem, we shall see that in this case our analytic approach 
is able to capture also the final beam profile edge details. 



B. Negative particles 

In the case of negative particles, the starting point is 
Eq. (pS)) (relevant under (strong) condition (01])), and 
again, it has to be examined with the object to trade 
the sum for an integral. First of all, it has to be minded 
that at n ~ 1 the hyperbolic arcsine arguments vary sig- 
nificantly, but at the same time they are large, whereas 
hyperbolic arcsine of a large argument is close to loga- 
rithm of a large argument: arsinhv ~ ln2t), and thus 
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FIG. 8: Mean volume reflection angle dependence on the crys- 
tal curvature, at = 400 GeV, for silicon crystal in orienta- 
tion (110). Points - experimental data from [l4j]. The line - 
prediction of Eq. (|72|) . with parameters evaluated as explained 
in the text. 



varies relatively little. On the other hand, in the do- 
main of large n the arguments of the arcsines vary lit- 
tle. Therefore, over the entire summation interval both 
sums involved may be approximated via integrals. Yet, 
first terms in the sums are singular functions of 
and therefore are better taken into account separately. 
Thereby, application of the Euler-Maclaurin formula to 
the first of the sums in Eq. (|46| gives 



arsinh 



n=0 



In 



25d {v^- 



d-26 



■In 



d-25 



V26diA^ 2 ^2Sd + 1) 

dn arsinh ^ 

!/(-)+! V2dd^yn 

1 - - 5 

— arsmh- 
2 



y^2M +ni„fi 



(73) 



1 



The next-to-leading order (derivative-related) correction 
term [13 to (1731) amounts 



arsinh 

12 dn 



2(5d(j/(-)+n) 



n— ninfi — 1 



24 (l-h 



(74) 

We will omit it because of the smallness of the numerical 
coefficient although, in principle, asymptotically it is 
also relevant (the same is true for all the higher deriva- 
tives, whose contributions enter with yet smaller coef- 
ficients (involving inverse factorial and Bernoulli num- 
bers)). 

Next, calculating the indefinite integral in Eq. ([75]) by 
parts. 



dn arsinh 



n ar; 



sinh^^ + a\/ d?- + n, (75) 
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we bring (j73l) to form 



"infl-l 

arsinh- 

n=0 ^ 



' + ninfl — 1 ) arsinh 



+ninfl-l) 



+- 



25d V 
1 + i^^"-* I arsinh 



+ vi-) + In 



d-25 



1, d-2d 
--In 
2 



— arsinh - 
2 



y^2(5d(l + iy(-)) 2 ^2^d(!/(-)+ninfi-l) 



(76) 



Now, in the thick-crystal limit riinfl oo, Eq. ((76 
simplifies to 



"■infl-l 

arsinh - 

n=0 



«inti»i ^ J25d 



.d- 2(5 



- (l-Fi^^-Mln 



d-25 



2(5d \ 



2Sd 



^25d{l + vi-)) 
d~2S 



2 (l + z.(-)) +ln , , 

(i-S) V25diA 



d-2S 



d- 2(5 



_ , d-25 /I , A, 

= ^ + ■ / 



d-25 
^/25dvi-) ' 



(77) 



Similarly, the second term in ()46p reduces to 

"infl-2 ^ + 5 

arsinh — ^ 
™=o y25d (i^(-) + n) 

, d + 25 /I , A, d + 25 

V25d V2 / ^25d(l + i.(-)) 



d-F25 



/ - -(l -f i'^"') -I- In - , 

25d 2^ ^ ^25dvi-) 



(78) 



Inserting ([77)1 and (1751) to Eqs. (pS)) and thereupon to 
(I52p . after a simple rearrangement one is left with the 




-l+'^-^ln(— In ) 



FIG. 9: Comparison of the exact formula (|52U46p for the nega- 
tively charged particle indicatrix (black solid curve) with the 
approximation ()79p (blue dashed curve) and approximation 
(|89p (red dotted curve) for R/Rc = 25. The central segment 
of the curves is strongly ^o-dependent. 



final result 



2Rc 
R 



■In 



In- 



e(l -//(-)) i^(-) 



1 



^(-) + ± In- 1_ 

2 J (l + i/(-)) 



(l-.(-))ln| 



(79) 



(the definition for ly^^^b) is Eq. (gOl)). Note that loga- 
rithmic asymptotics of this expression at i'^^^ — > and 
at — >■ 1 agrees with the general law ()55dP . 

The exact ([52]) vs. approximate ((79|) expressions for 



the indicatrix 



(6)) are compared in Fig. [9l 



1. Differential cross-section 

To deduce the observable differential cross-section from 
the available indicatrix, we have again to issue from 
Eq. (j65p . To some extent, the same procedure as for 
positively charged particles applies here, leading to rep- 
resentation 



dX 



dO^ 



2 

E 



dO^ 



1 



\duY^/dh\ 



(80) 

where ('^^ ■*) is the solution of Eq. (|i(It . X]j=i ac- 
counts for the existence of two roots for equation 6'v.r. — 
dv.r. (i^*' with the function ([79| . Again, upon aver- 
aging over ^0 under condition ([69]) (cf. Appendix), one 
obtains: 



1 

d \de. 



dX 



E 



di^ 



(-) 



de. 



i. e., averaging over b (and a tiny interval of ^o) reduces to 
averaging over (transverse energy). Note that when 
becomes a uniformly distributed random quantity in 
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a unit variation interval, for each given b, and thus ^v.r., 
one can always unambiguously tell whether 
or u^~^ = v\ \ so the probability normalization is con- 
served under conditions of summation over the branches. 

A technical distinction of the case with negative parti- 
cles is that equation ([7^ can not be resolved with respect 
to v^~^ in an explicit and exact form. Of course, it can 
be easily solved numerically; the differential cross-section 
so evaluated is shown in Fig.[TUl by a solid curve. On the 
other hand, it is also useful to pursue an analytic but ap- 
proximate approach, based on different approximations 
in different regions of J^^^-*, to which we yet pay some 
labor. 



6»,. dA 




2. Asymptotic evaluation of the final beam profile 

In Eq. (17^ at typical z^'-"-' the leading term is the last 
one, where j/^"' is multiplied by a large logarithm. Be- 
sides that, in the domain of i^^"-' close to 1 the first log- 
arithm in (j79p becomes large, too, and a minimum of 
function 0v.r. (^^^"•') develops, corresponding to onset of 
a rainbow scattering. As for the second logarithm in ([7^. 
which raises at i'^~'> — >■ 0, it does not lead to formation of 
a dependence ^v.r. ('^'' minimum - on the contrary, it 
makes the dependence steeper, and in the area of its sig- 
nificance the differential cross-section is small (exponen- 
tially). So, for the differential cross-section description it 
basically suffices to consider only two regions: the region 
where the last term of ([79]) dominates, and the region 
where the last term and ^ In ^(~)) ^^^P compet- 
ing. On the 0v.r. axis the mentioned regions are adjacent, 
and conjointly they should give almost the full picture of 
the differential cross-section variation. For completeness, 
one may consider also a third, asymptotic region of the 
the differential cross-section tail (orbiting region) , where 
the first and the second logarithms of (17^ dominate. 

a. Rainbow region. The value of v^"^ corresponding 
to the rainbow angle is to be determined from condition 



90v.r. 



= 0, 



(82) 



which in application to expression ([75]) gives 



FIG. 10: The 0o-averaged differential cross-section of nega- 
tively charged particle scattering, at R/Rc = 25 (solid curve). 
The axes scales are chosen so that areas under the curve is 
unity, as the total probability. The left-hand dashed curve is 
evaluated by the explicit approximate formula (|86p . the right- 
hand dashed curve - by the explicit approximate formula (|93|) . 
It is observed that those approximations actually overlap. 



tion ^v.r. (^^' ') up to a quadratic term 
R 



2Rr 



1, /16, e^/^i? 
In ^ m 

^ 2Rr 
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In' 



2Rc 
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Now, expressing the pair of roots Vj (6'v.r.) from the 
quadratic equation ([55)) . one derives by formula (|M|) the 
behavior of the 6'o-averaged differential cross-section in 
vicinity of the rainbow angle: 



Br I dX 



d \d9v.i 



VR/2Rc 



(86) 

The domain of applicability of this approximation is de- 
termined from Eq. (|85p by demanding the third-order 
term to be small compared with the second-order one: 



,(-) 



1 - v^-^^ ly^-^ 



In 1 



The approximate solution of Eq. 



« 1- 



1 



2 In 



2R^ 



O 



R 



= ln— »1. (83) 



IS 



1 



In' 



3 e-V-ifl 
2R^ . 



In vicinity of the found point ^ we may expand tunc 



(84) 



1 1 — ml — m — — I ^ 1 (rambow region). 

(87) 

Function (jSG]) is shown by the left-hand dashed curve 
in Fig. [lOl 

6. Typical angle region. On the other hand, if z^*^"^ 
is not too close to 1, i. e. under condition 



R , /16, e5/4i?\ 



- 1 1 In ( — In ) > 1 (typical angles) 

"■^ 2Rc J 
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opposite to ((87|) . all the terms in (|79p except the second 
logarithm (singular at — >■ 0) may reasonably be ap- 
proximated by their Taylor expansions up to linear terms 
- say, about the unit interval midpoint v'^-^^ = \ [soj : 



. 2RJ ViR 11 (_), 2i? 

1 In^ 'a\ r-rr 



(89) 

Then, for determination of the inverse function 
v'^~^ (Ov.r.) one obtains the following simplified equation: 



1 



R 

2Rc 



11 



In 



In 



2R 



SeV^Rr 



Here the r. h. s. is a monotonic function of 
so, in contrast to the exact equation (j79p . the simpli- 
fied equation has only one root (the second, lost 
root gives a small contribution to the differential cross- 
section). The solution to Eq. ([TO]) expresses through the 
Lambert (product log) function W{s) defined as a solu- 
tion to equation In s = In -I- W, and incorporated in 
many computational software packages: 



1 



In 



^/y3i?ln3^ 



2R 



Rc 



(91) 

(Mind that 9y,^ J6c is negative and close to —1). From 
the definition of W , its asymptotic behavior in different 
regions derives as 



In, 



S < 1 
S > 1, 



(92) 



and its derivative 

W\s) = 



s{l + l/Wis))' 



Therefore, over the typical angle region the 6'o-averaged 
differential cross-section is cast as 



dX 



d \de^ 



di^^- 



de^ 



R 



1 



2i?cln3^(l 



i/wy 

(93) 



the argument of W being the same as in Eq. (|9ip . 

The full width of the differential cross-section as can 
be inferred from Eq. ()93p . and is obvious already from 
Eq. ((7^ . amounts 

2R R 

A^v.r. ^ 0c In — (negatively charged particles) , 
R Rc 

(94) 

which at i?/i?c > e'^ w 23 is larger than width (|64p for 
positively charged particles. 

The behavior of function (|93| in region (|88l) is shown 
in Fig. [TU] by the right-hand dashed curve. It describes 
the actual distribution quite accurately. In the orbiting 
region asymptotics, however, Eq. (j93|) errs by a factor of 



e^/e/VS « 1.3 (see Eq. ^ below), so apphcability of 
(|93l) is restricted by the bound 



R 

2Rr 



1 



^ Rc 



(typical angles) 



in addition to 

One may anticipate formation of a shoulder in the dis- 
tribution ((93)) at sufficiently large R/Rc, since then the 
argument of W can be large and 1 + l/W ^ 1, making 
the differential cross-section 6'v.r. -independent. However, 
W{s) reaches 5 (» 1) only at s ~ 700, so a shoulder in 
{dX/d9y,r.)gg may form up only at R/Rc > 100. 

c. Orbiting region. Finally, asymptotically large 
6'v.r. are generated in the regions — >■ and — 1, 
where there are logarithmically rising terms in the rela- 
tion ((79)) . Examine first the region 0. In this 
limit one may let = everywhere in ()79p except in 
In ijzrf, and the reduced equation 



1 



is easily solved: 



eRn 



2Rc R 



(95) 



Differentiating that expression with respect to (?v.r., the 
differential cross-section asymptotics results as 



6r I d\ 



d \d9v.i. 
holding if 



i?2 



2ei?2 



2Ra{-^ ; (orbiting region), 

(96) 



R f I- \ , R , , . . .X 
— ( + 1 1 > In — (orbitmg region). 

As for the contribution from region z^^^-' — )■ 1, which 
can be treated in a completely analogous manner, it 
equals 

e^Rc 

and is definitely negligible compared to ((96)) . so Eq. ((96)) 
is the complete asymptotic result. 



3. Comparison with experiment 

The only measurement of negatively charged particle 
volume reflection is jj]. It gave, for Si (110) orienta- 
tion, R/Rc ~ 70, the volume reflection angle mean value 
|(6'v.r.)| — 0.666'c. This is significantly smaller than our 
expectation (and other simulations, as quoted in (isj ) 
|(^v.r.)| ~ Sc (yet minus ^ 8% correction for finite Rc/R, 
which does not matter anyway). Of course, if we treat 
I-Finaxi as an adjustable parameter of the model, we might 
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achieve agreement with the experiment, but a physical 
recipe for this is yet lacking. 

As for numerical calculation results for the final beam 
profile, there is lack thereof for negatively charged par- 
ticles at R ^ Rc and free of multiple scattering (and 
with boundary conditions congruent with typical exper- 
imental ones) [3l|. We hope to see such results in near 
future. 



V. HIGH-ENERGY PASSAGE LIMIT 
(PERTURBATIVE DEFLECTION) 

In conclusion, we will briefly comment on behavior of 
the function 0v.r.(&) in the opposite, high-energy limit. 



d R 



(97) 



Thereunder, the deflection becomes perturbative (and 
better viewed in Cartesian coordinates, without the ref- 
erence to a centrifugal force notion), and for positive and 
negative particles it must be equal in magnitude but op- 
posite in sign. That is conflrmed by Figs. |3] and Fig. [TTl 

The speciflc expression for the dependence 0v.r.{b) 
this limit was obtained in [l5| (Eqs. (18-19)), in the 
Cartesian coordinate framework; 



9.A0o,b) ^ 0Born{Oo,b) 



1 

E 



dzF{b, z) (98a) 



^^V2Rdj 1 r 1 



Rn2 , b 



(98b) 



(again, the braces indicate taking the fractional part). 
It involves the Hurwitz (generalized Riemann) zeta- 
function at a negative value of its flrst argument, which 
may be defined, e. g., as a contour integral [13] 



r(i -g) 



(0+) gQ-lgfS 



1 



ds 



(99) 



along a Hankel path (S^j- Function (j98b[) is shown in 
Fig. mby dotted line. Note the identity C(-5,0) ^ 
i,l) (= C(— ^)): whereby function (j98bl) is every- 
where continuous, but its derivative breaks at point 
6 = 1 _ /Ml 

d 2 \ 2d Jj- 

The second argument of C-function in (j98b|) allows for 
a physical interpretation |16| : 



(100) 



is the particle impact parameter at a depth where the 
particle straight trajectory becomes tangential to the 
bent atomic planes (actually, trofl)- Should one pass to 
6rcfl instead of b, the dependence on 9o disappears com- 
pletely: 







/' 










0.5 


\ /■ 


1 






-0.5 






4f(-l/2) 





FIG. 11: Comparison of impact parameter dependencies of 
the reflection angle for positively charged (solid curve) and 
negatively charged (dashed curve) particles, at R/Rc = 1/20. 
Dotted curve - approximation (|98bp . for positive particles. 
The functions displayed may be continuously periodically ex- 
tended beyond the interval ( — |, ^)- The position of the frac- 
ture at ^ = i — I I corresponds to tangency of the parti- 
cle near-straight trajectory to one of the bent atomic planes. 



This contrasts with the case R^ Rc investigated in the 
previous section, where the 0o-dependence yet remained 
in a specific, casually located rainbow peak. 

Here we will not contemplate demonstrating from ex- 
pressions peiSObI) that the limit of 6'v.r. is indeed ((98b| . 
We only note that to this end one must implement signif- 
icant cancellations between the pair of terms under the 
sum sign. One also notes that the second argument of 
(■-function in (j98b|) is just the limit of i/(+): 



(0O,&)=^v.r. (&rcfl). 



(101) 



For negative particles the route to (j98bp is a bit more 
intricate. In Fig.[TT]we numerically compare approxima- 
tion (I98bp with the exact result - the agreement is quite 
convincing. 

The comparison with the perturbative scattering pat- 
tern may also give an insight into the origin of the volume 
reflection phenomenon. The average of function (j98|) over 
the impact parameter b turns out to be strictly zero (see 
[l6| ) - hence, in the high-energy limit signatures of the 
volume reflection completely disappear. That owes to the 
fact that J'^^f^2 dbF{b, z) = at any given z (regardless of 
the crystal bend). In contrast, at R> Rc the trajectory 
may not be viewed as merely straight [ssj . and the par- 
ticle distribution is non-uniform over the crystal volume. 
In fact, "shadowed regions" not filled by the particles 
may appear at the inner side of bent potential ridges, in 
which the force acts in the positive direction. The deficit 
of a positively directed force on the beam then leads to 
the negative sign of the particle beam mean deflection 
angle. 



17 



VI. EXPERIMENTAL OBJECTIVES AND 
PARAMETER OPTIMIZATION 

The alleged application of volume reflection is for high- 
energy particle beam extraction from accelerator beam- 
lines, when that is not manageable locally with labora- 
tory magnets (intra-crystal forces are much stronger). 
Thereat, the beam parameters (energy and angular di- 
vergence) are definite, while the crystal parameters have 
to be optimized in order to attain a suitable deflection 
quality. 

Control of beam divergence accompanying the 
deflection. The mean deflection angle 



= Y ^ ~ 2lAirad 
(for Si (110)) 



100 GeV 



E 



(103) 



for a silicon crystal of orientation (110) depends only on 
the particle energy. The next parameter to care about is 
the angular spread acquired by the beam at the exit from 
the crystal. Neglecting the incoherent multiple scatter- 
ing (to be estimated below), the spread A6'v.r., given by 
equations ((64|) and ((94|l . for a given beam energy and 
crystal material depends only on the crystal bending ra- 
dius. To derive a criterion for beam complete deflection, 
one may demand that the bulk of the dispersed beam 
(its both "edges") be deflected to the same side. That 
implies: 



R > ARc (for positively charged particles) (104) 



and 



R 



R > 2Rc In — . (for negatively charged particles) 
Rc 

(105) 

Relation (|104p with the factor 4 was found previously by 
Maisheev Q based on numerical simulation studies for 
protons. Our paper, thereby, offers a formal justification 
for that empirical relation, although within a framework 
of a simplified model. The emerging ratio 



q = 



max|6lv.r.| j iS; for pos. charged particles 



A(?v.r. 



^ R for neg. charged particles 

(106) 

quantifies the steering quality. But for the deflection to 
be neat, one should have a substantial q 

q>3 (neat deflection), (107) 

entailing for the crystal curvature radius 

R>l0^l5Rc. (108) 

In demonstrational experiments Q the latter require- 
ment was marginally satisfied. 



Sufficient crystal thickness. The crystal thickness 
(i), it to be kept as low as possible in order to minimize 
the multiple scattering. However, a lower bound for L 
results from the requirement that the volume reflection 
has space to develop, i. e., the "thick crystal limit" holds, 
in the sense of Sec. IIII CI Thereto, the crystal bending 
half-angle ^ needs be larger than the critical angle 9c, 
which is L-independent. Thus, 



L > 2R9c 



R 

Rc 



E 

Wo 



J_3.6Mm^^ 



(109) 



(thick crystal, volume reflection saturation) . 



If we regard here = q as fixed, in general L grows 
with the energy At E - lO^-^-lO^ GeV (RHIC, Tevatron) 
the minimal thickness amounts only to (0.036-^0.1 mm)q. 
Note that crystals as thin as 30 mm are manufacturable 
(as described in [H). For E - lOTeV (LHC) and q - 3 
minimal L reaches the value of 1 mm. 



Reduction of multiple scattering. To quantify 
the impact of incoherent, random multiple scattering, we 
have to evaluate the characteristic ratio [3J| 



gmult 



13.6 MeV 

^/2E 

0.47 



93.6 mm 



E 

Wo 



100 GeV L 



E 



(Si (110)). (110) 



If the latter square root does not exceed unity (in ex- 
periments Q it is ~ 1), the multiple scattering does not 
spoil coherent beam deflection. We conclude that for 
-E > 100 GeV crystal lengths up to 2 mm are multiple 
scattering safe, i. e., ^^""/^c < 1- 

Besides aggregate deflection, it would be interesting to 
experimentally investigate the intrinsic volume-reflected 
beam shape, and in particular to check the shape depen- 
dence on the particle charge sign (cf. Figs. [71 and [TU]) . 
The main problem here is that at q ^ 1 the final beam 
half divergence ^A^v.r. is 2q times smaller than the mean 
deflection angle, and so is sooner overtaken by the multi- 
ple scattering, making the profile Gaussian and particle 
charge sign independent. For this not to happen, one 
needs condition 



amult 

2q^ « 1, 



(111) 



with 6'j;'"'V'^^v.r. to be inferred from piU|) . Hence, for 
the present purpose we should not strive for large q, 
granted that the final beam profile is not very sensitive to 
q at q > 1. So, q ~ 1.2 seem to be good enough. Equally 
well, in order to raise the angular resolution we should use 
moderate energies. Say, E = GeV is ultra-relativistic 
enough. To reduce multiple scattering, we can take a thin 
crystal with L = 30 /^m, which marginally satisfies (|109l) . 



This gives 



3 • 10-^ rad, M,, 



2.5 • 10"^ rad, and 
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— ~ 0.3 (small enough). But one has to control ini- 
tial particle impact angles with an accuracy a few times 
better than A9y,r.- This may be difficult to achieve via 
initial beam collimation alone, so one may need to apply 
event selection procedures (cf. Q). 



VII. SUMMARY 

Based on the model of a purely parabolic continu- 
ous potential in a bent crystal, we have gained a lot 
of information about the volume reflection phenomenon, 
for cases of positively and negatively charged particles. 
First, we have obtained an explicit expression ([^51) for 
particle trajectories. From the solution for trajectory, in 
particular, we have derived the particle final deflection 
angle as a function of the particle impact parameter and 
energy, in form of sums (j361 150b|) . Asymptotic behavior 
of those sums at R ^ Rc was explored, and asymptotic 
values for the volume reflection angle were found. They 
equal: — f for positive particles, and —9c for negative 
particles. This agrees within ~ 20% with the existing 
results of numerical simulation using more realistic con- 
tinuous potentials 0, and with experiment for positive 
particles Q (though there is an indication of worse agree- 
ment for negative particles Ql). 20% is about the same 
accuracy as for approximating the continuous potential 
by a parabola. Yet we have evaluated the next-to-leading 
order correction in parameter Rc/ R, which depends on 
the impact parameter, and, by averaging over impact pa- 
rameters, we determined asymptotic shape of the final 
beam. This in particular yields the mean volume reflec- 
tion angle dependence on Rc/R, which appears to be 
linear - in general agreement with experiment (l^ | (see 
Fig.il). 

In course of investigation of the final beam shape, we 
have discovered various singularities in its profile, which 
moreover are particle charge dependent. First of all, we 
had to deal with the problem that, in principle, the final 
beam profile may contain a visible admixture of boundary 
dependent effects ("randomly" located peaks). However, 
we have proved the statement that boundary effects get 
completely erased in the differential cross-section aver- 
aged over a tiny interval of incident angles Oq (condition 
([691) '). 01"' analogously, due to a bit of multiple scatter- 
ing before the volume refiection region. Therewith, the 
averaging over impact parameters becomes equivalent to 
averaging over parameters z/^*^ (i. e., transverse energy), 
and we were able to analytically deduce the final beam 
profile for positive and for negative particles. For neg- 
atively charged particles it is asymmetric, exhibiting a 
spike on its outer edge, corresponding to the rainbow 
scattering, and an exponential tail on the inner side, cor- 
responding to orbiting (Fig. [TUl). For positive particles, 
the final beam has a rectangular profile (Fig. But in 
actual practice, with the account of continuous potential 
smearing in vicinity of the atomic planes, one expects 
appearance of a weak rainbow spike and orbiting tail for 



positive particles, as well. 

Towards practical applications and further experimen- 
tal investigations, we have made a few numerical esti- 
mates. They indicate that for usage of a bent crystal as 
a coherent beam deflector, one needs a relation between 
the main parameters 



< 



E 



< 



R /20^rad 



1mm 100 GeV m' V ctq 



(ctq is the r.m.s. angular deviation in the initial beam). 
The better those inequalities are met, the higher is the 
deflection quality. If one becomes interested in investiga- 
tion of the flnal beam intrinsic shape, generated by the 
continuous potential alone, those inequality must be sat- 
isfied strongly, but minding existence of technical lower 
limits for L and cto. This suggests an optimal energy 
about 50 GeV; experiments are to be carried out simul- 
taneously with particles of both charge signs (e*, tt*). 

There are many respects in which the model solution 
described herein can be improved. First of all, it is 
straightforward to add to the simple parabolic poten- 
tial a second parabolic section - either to round off the 
potential in vicinity of atomic planes, or to describe the 
potential of Si crystal in (111) planar orientation, which 
is of practical importance, too. As a next step - at least 
a perturbative account of incoherent scattering processes 
is desirable. But at the same time, even in the present 
form, the theory (trajectories derived in Sec. seems 
suitable, e. g., for study of electromagnetic radiation 
emitted by a volume-refiected particle. 



Appendix A: Formal procedure of ^o-averaging 



In Sec. IIV Al we had obtained, for positively charged 
particles, the scattering differential cross-section in form 



dX 



R 



de. 



2tt0cR, 



E 



1 



(Al) 

Our objective now is to straightforwardly compute the 
sum involved hereat for the specific function i'^~^\b,9a) 
given by Eq. psp. first for an arbitrary and then 
average it over 9o, in order to justify our assertion that 
combined averaging over b and (within tiny interval 
(|69p ) is equivalent to averaging over 

To begin with, let us find the roots bm explicitly. 
Eq. (|55l) is equivalent to 



.2/12 



+ {b + Sr-{l-6)' ^ 
2Sd 



(A2) 



with m an integer, solution of which is straightforward: 



b^{v^+\9o) = -5±\j25d +m)+(--Sj - T^9l 

(A3) 
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Here the sequence of m begins with a smaUest integer mo 
at which the radicand in (|A3p is yet positive, viz. 



mo = 



25d 



1. 



(A4) 



The upper hmit of m in Eq. (IA3p equals to the largest 
integer at which yet & < f (for a branch with the "+" 
sign in front of the root in Eq. (jA3p ). and 6 > — | (for a 
branch with the "— " sign). That yields, correspondingly, 
values 



^inax 1 



Llo _ ^(+) 



(A5) 



26d 

'71inax2 = rn^axi + l- (A6) 
The number of terms in the sums from (IA4I) to (IA5[ IA6I) 



is large: 



TOmaxl,2 -mO=^+ 0{1) = — + 0(1) > 1- 



86 



4i?c 



Next, values of derivative in points 5m are easily 

evaluated, noticing that the fractional part operator in 
i/*^"*"' is inconsequential for derivatives |35| . Differentiat- 
ing (|A2p gives 



d6 



5d 



Identity 



db 



dm 



(A7) 



(A8) 



suggests that upon substitution to (|Aip one may expect 

1 ^^'T'max 1 



E + E 

^,m— mo m— mo y 



2 E 

m— mo 
fd/2 



db_ 

dm 



2 db = d, (A9) 

"'0 



as we anticipated in Sec. IIV A[ but to accommodate the 
dependence on ^oj need to carry out the calculation 
more precisely. 

Through (IA3[ lAll) . the differential cross-section as- 
sumes the form 
dX{d ) R 



dOy.,. 2T:ecRc 



Sd 



E + E , — 

m=o rn=o J ^2Sd + m) + {'I- S) -T^Ol 



2tt6 
with 



R/4R. 
It \ ^ 

it ^ 



'v.r. J, fo 



'o y^m-|-a(i^(+)(6'v.r.),6'o) 



(AlO) 



= 1 - 



26d 

{b{iy^+\e^.,.),eo) +6y 

25d 



.(All) 



Outside the interval 



2 



^ > (A12) 



there are no roots to equation i'*-^-' (6'v.r.), so the differen- 
tial cross-section vanishes as an empty sum. 

At a large upper limit the sum in (jA10[) grows as 
^R/Rc, whereas the difference 



F — ^ 

^0 ^"^ + " 



R 

Rc R/Ra 



-, a 



(A13) 



tends to a finite limit (it may be categorized as the Hur- 
witz, or generalized Riemann, zeta-function, see also def- 
inition ((M)) ). Thereby, we may cast (|A10p as 



dXi0. 



v.r. , 170 J 



de^ 



d R 

2^r~Rr 



2' 



k(ev.r.,0o) 



(A14) 

Here ^jRc/R ^ 1, and the unity in parentheses (|A14p 
corresponds to the result anticipated in (|A9p . But one 
should take care that at a — function C blows 
up as So, the correction in the parentheses in (IA14P 
can not be regarded as everywhere small. 

Now, we turn to the issue that a is ^o-dependent. 
When ^0 varies (at 0v.r. fixed) even in a narrow interval 
(see l[67ll). a uniformly and repeatedly scans its defini- 
tion interval from to 1. Hence, the averaging over ^o 
is equivalent to the integration over a from to 1. By 
virtue of the property 



daC, ( 7:, a 1 = 



(A15) 



checkable from definition (IA13p . we have 

d\ \ f\ dX 

da- 



d9v.r. I Jo d0v.r. 

Rd 



4Rr 



Thereby, we arrive at equation (I7ip . without any correc- 
tions ^ Rc/R, q.e.d. 

For negatively charged particles there arises a sum sim- 
ilar to that of (|Aip . though the dependence v^^^{b) dif- 
fers a bit from v^^\b). Nonetheless, the averaging proce- 
dure is completely analogous, the non-averaged differen- 
tial cross-section equals ((5T|) times the parentheses factor 
of (jA14p . Upon the ^o-averaging, through (|A14p again, 
we arrive at Eq. ((5T|) . 
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and thus coherent action of atomic forces over long dis- 
tances. Still, quantum effects may be viable in a special 
case when the particle transverse energy is very close to 
the height of an atomic potential barrier (the author is 
indebted to A.V. Shchagin for pointing this out in private 
conversation), but in any case, the classical calculation 
has to pave the way. 

[25] Equation U]) in itself may apply to a non-ultra-relativistic 
motion, too, provided E includes the particle rest energy 
(so, in the non-relativistic limit E — >■ m). But physically, 
in order not to complicate the analysis, we confine our- 
selves to (the most important) ultra-relativistic case in 
the present paper. 

[26] Strictly speaking, the trajectory will not be exactly sym- 
metric with respect to t^ca because the distances from t^cfi 
to the crystal boundaries are in general different. How- 
ever, contributions to Srefl from crystal regions away from 
trcfl are supposed to decrease sufficiently rapidly, and one 
expects existence of a "thick crystal limit" of ^rcfi, rele- 
vant in actual practice - see Sec. IIVI In this paper, we 
content ourselves only to the "volume" contribution ((SSj 
to Srcfl and do not study any boundary effects. As we shall 
see later (Sec. llVl and Appendix), however, the omission 
of boundary effects requires certain care. 

[27] The verbal description of the given effect for negative par- 
ticle case is contained in 0] (end of Sec. 3), and a graph- 
ical illustration thereof appears in (Fig. 5, trajectory 
2). But relation with the general concept of orbiting (spi- 
ral scattering), as formulated in [l^li was only noted in 
[13 ]. We don't go into discussion of possible experimental 
significance of orbiting here. 

[28] Equation (|55a|) agrees with our definition of the 9 angle 
and differs in sign from the conventionally defined angle 
in a centrally-symmetric field. 

[29] In this subsection the sign alternative in f ±5 corresponds 
to dealing with the first or with the second of the sums 
in (f36bl) . 

[30] More elegant expressions, at the same time accommodat- 
ing for a correct asymptotics in the orbiting region, result 
when the Taylor expansion is carried out about the point 
— 0, but unfortunately, that approximation would 
be definitely less accurate at typical 6v.r.. 

[31] Paper does not exhibit negative particle results for 
-R > 7?c. In the pioneering work Q] the boundary condi- 
tions of are different from ours and from typical experi- 
mental arrangements: the authors consider particles en- 
tering the crystal face parallel to the bent atomic planes 
(a "half volume reflection). That should produce a dif- 
ferent final beam shape (with no rainbow singularity and 
more particles in the "orbiting" tail), but nonetheless. 
Fig. 4b of [3] is visually rather similar to our Fig. [TOl 

[32] Conventionally, the Hankel path is defined to begin in the 
complex s-plane at — oo (args = — vr), encircle the origin 
in the positive direction and return to — oo (arg s — -\--k). 
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[33] In the literature, sometimes, one meets an interpretation 
of the volume reflection phenomenon right in terms of 
a straight trajectory tangential to the bent crystalline 
planes in some point. Such an interpretation, although 
satisfactory for a symmetric installation of crystals with 
a large bending angle (~ 29o ^ 9c), may be misleading 
for understanding of the underlying particle dynamics. 

[34] For the (rms, plane) multiple scattering angle upon the 
particle over-barrier passage we crudely apply a formula 



for the scattering angle in an amorphous target made of 
the same material (silicon): eij^"" = "^"^ 
(as quoted in [l9j). 
[35] Strictly speaking, differentiation of finite discontinuities 
will give (5-functional terms, but they will be impercepti- 
ble when inserted to the denominator of equation (|65[). 



